Nested Quantum Error Correction Codes 
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The theory of quantum error correction was established more than a decade ago as the primary 
tool for fighting decoherence in quantum information processing Although great progress 

has already been made in this field, limited methods are available in constructing new quantum error 
correction codes from old codes. Here we exhibit a simple and general method to construct new 
quantum error correction codes by nesting certain quantum codes together. The problem of finding 
long quantum error correction codes is reduced to that of searching several short length quantum 
codes with certain properties. Our method works for all length and all distance codes, and is quite 
efficient to construct optimal or near optimal codes. Two main known methods in constructing new 
codes from old codes in quantum error-correction theory, the concatenating and pasting, can be 
understood in the framework of nested quantum error correction codes. 

PACS numbers: 03.67.Pp, 03.67.-a, 89.70.Kn, 03.67.Lx 
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Quantum computers offer a means of solving certain 
problems, including prime factorization, exponentially 
faster than classical computers One of the main 

difhculties in building the quantum computers is that 
the quantum information are fragile to decoherence. To 
solve this problem, the quantum error correction codes 
(QECC) are designed which provide us an active way 
of protecting our precious quantum data from quantum 
noises. Since the initial discovery of quantum error- 
correcting coeds (QECC), researchers have made great 
progress in codes construction. But for large number of 
qubits there has been less progress, and only a few gen- 
eral code constructions are know. The main difficulty 
in constructing long codes is that we need to choose the 
suitable QECC from a large quantity of quantum states 
whose number in general grows exponentially with the 
length. This random search method will be intractable 
for long codes. However, long QECC is indeed necessary 
when we can control more and more qubits to realize 
the scalable quantum computation. If we can divide the 
searching of QECC into several easier steps, then the dif- 
ficulty in constructing long QECC will be reduced dra- 
matically. In this work, we propose a general and easy 
method to construct the nested QECC. 
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I. CODES CONSTRUCTION FOR DISTANCE 3 

A. An optimal [10,4,3] QECC 

Our main idea to construct long QECC from known 
short length QECC is to ensure that all syndromes for 
one-qubit error of the long QECC are different. For ex- 
ample, we may consider to repeat a [5, 1, 3] QECC twice, 
automatically, the syndromes are also repeated twice. 
Since the syndromes for [5,1,3] are different, the next 
step to construct a long QECC is to find additional gen- 
erators so that the syndromes in each block can be distin- 
guished. So we consider to put a code of length 2 whose 
syndromes are also different under any block to achieve 
it. We can find the additional generators that take the 
following form. 



1 
1 1 



1 1 
1 



(1) 



We name it subcode {2,2}, i.e., the length is two and 
the code has two generators. To ensure that it is a valid 
QECC, we also need to confirm that each generator of the 
above matrix, (YZ) and (XY), which will be repeated 5 
times in the new code commute with each other and also 
commute with the original generators of [5, 1, 3] which is 
named blockcode. Suppose in GF{2) each generator of 
subcode is written as {a\b) = (aia2|&i&2) and for block- 
code as (c|d) = (ciC2C3C4C5|did2rf3c'4c'5), then commuting 
condition for this construction will be 



(ai -I- a2)(di + d2 + d3 + d4 + d^) 

+ {bl + fc2)(ci + C2 + C3 + C4 + cs) = 



(2) 



(see (|A19[) and (jA20[) for details about how to obtain 
this equation). However, we may consider a stronger 
constraint condition, for example, (ci-|-C2 + C3-|-C4-f C5) = 



2 



and (c?i + ^2 + + 1^4 + rfs) =0 which is for the 
blockcode. Really we can find a [5, 1, 3] code that satisfies 
this condition and the code takes the form 



/ 1 1 

10 1 

10 10 

Vo 1 1 



110 

110 

1 1 

1 1 



(3) 



So stabilizer codes constructed from [5,1,3] and {2,2} 
presented in the following constitute a new QECC. 



(4) 
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This is a [10, 4, 3] code and we know it is optimal. 

B. Other optimal QECC constructed from [5,1,3] 
We can also find subcodes {3, 2} and {4, 3} 

(5) 



(6) 
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10 1 
11 
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We consider the same stronger constraint condition for 
blockcode. From [5,1,3] code, the new QECCs con- 
structed are [15, 9, 3] and [20, 13, 3] which are all optimal. 
Their explicit forms are presented in the appendix (see 
(EStI) and jJMI). 



C. QECC of length 25: nesting two [5,1,3] codes 

A straightforward method for a length 25 codes from 
[5,1,3] code is the concatenated QECC [J- Here our 
method is quite different from the concatenation. To 
construct a QECC for n = 25 from code [5, 1,3], a sim- 
ple method is to let both the block-code and the subcode 
be a [5,1,3] code. Then commuting condition for this 
construction will be 

(ai + a2 + as + 04 + 05) (di + ^2 + c?3 + c?4 + d^) 

+ (61 62 + 63 + 64 + 65)(C1 + C2 + C3 + C4 + C5) = 0. 

(7) 

However, we may consider a stronger constraint condition 
as (01-1-02-1-03-1-04-1-05) = Gand {C1+C2+C3+C4+C5) = 0. 
Then we can find a [5,1,3] code for both blockcode and 



subcode that satisfies this condition and the code takes 
the form 



1 1 

10 1 

10 1 

1 1 



111 
10 10 
1110 1 
1110 



(8) 



The newly constructed QECC is a [25, 17, 3] code, see 
(|A9|) which is not optimal. To construct the optimal code 
for 25-qubit which should be a [25, 18, 3] code, we need a 
{5, 3} code to be the subcode. Still we use the stronger 
constraint condition {ci + C2 + C3 + C4 + c^) — and 
{di + d2 + ds + d4 + d^) = for the blockcode, then there 
will be no extra conditions for the subcode, so any {5, 3} 
code can be used. One example of an optimal [25, 18, 3] 
code constructed is presented in the appendix ,see (|A29|1 . 

We may notice that that syndrome (0000) of code 
[5, 1, 3] which is not used in constructing code [25, 17, 3]. 
With this property in mind, we can repeat this syndrome 
5 times and then nest a [5,1,3] as the subcode to this 
syndrome. Then a code [30, 22, 3], see (jA12p can be con- 
structed, this is a near optimal QECC. Continuously we 
can construct a [35, 27, 3] code, see (jAl4[) which is also 
near optimal. 



D. Method to constructing perfect QECC 

Here we present our nested method in constructing one 
class of codes, perfect codes |4|, which are optimal since 
it achieves the Hamming bound. 

First we need a class of sub-codes, we name them as 
raw perfect-constructing codes, which are one kind of 
{2'^ - 1, k} QECC whose all (2'^ - 1) one-qubit errors 
take all the 2^^ — 1 syndromes that the code can provide 
except for syndrome (00 • ■ • 0), and so do the cr^ and ay 
errors. Suppose any generator of a {2'^ — 1, fc} raw perfect- 
constructing code is (ci,-- - , C2fc_i|(ii, • • ■ ,c?2fc-i), one 
feature of those codes are (ci + • • • + C2fc_i) =0 and 
(di -I- • • • -I- ^2^-1) = 0. So if we choose a raw perfect- 
constructing code as subcode, all commuting conditions 
satisfied automatically. So there will be no extra condi- 
tions for the blockcode which can make the code search- 
ing much easier. 

We have the following results: 

(1), A [2^^, 2*^ — fc — 2, 3] code can be constructed by nest- 



ing a "code" 



as blockcode and a [2*^ — 1, 2'' 



1 — fc, < 3] raw perfect-constructing code with syndrome 
(00 • • • 0) as subcode together. 

(2), A [(2'=+'=' -l)/3,(2'=+'='-l)/3-(fc-Hfc'),3] perfect 
code can be constructed by nesting a [(2*^ — l)/3, (2'"' — 



l)/3 - fc, 3] perfect code, a [2*= -1,2'= - 1 - fc', < 3] raw 
perfect-constructing code and a [(2*= — l)/3, (2^= — l)/3 — 
fc',3] perfect code together. 

We take the constructing of [16, 10, 3] code perfect 
code as examples and present the detail in appendix, see 
(1X36)) . 
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Gottesman's stabilizer pasting of distance 3 code is also 
easy to be understood in our theory, but our method is 
more general for two reasons: (i) blockcode of [712, S2] can 



be any code not only 



which has been shown in 



constructing [30, 22, 3] code and [35, 27, 3] code and per- 
fect code presented in above; (ii) the syndromes of block- 
code of [ni, si] need not to be 00 • • • only. In Ref . pH . 
Yu et al. present a construction of optimal [37, 30, 3] code 
which is [37] = [2^] > [5], here we give another construc- 
tion which is more powerful, the detail can be found in 
appendix, see (|A45ll . 



II. NESTING CODES FOR ALL DISTANCE 

Not only for distance 3, our method works for all dis- 
tance QECC. Suppose there are two copies of a subcode 
whose generator matrix is [A], the easiest way to distin- 
guish their syndromes is to connect them in block diag- 



onal matrix as 



0^ since the logical operators of 

this construction are only constructed by logical opera- 
tors of these two codes, and the hardest way is to con- 
nect them parallel which is [A A). For the hardest 



way, the problem of constructing distance 3 codes have 
been solved above, and problem of constructing larger 
distance codes is still open. But for the easiest way, codes 
construction for all distance can be solved which will be 
presented in below. 

Suppose there are m subcodes, i.e., [ni, /ci, di], 
[n2,fc2,(i2], [rim, km, dm], that are con- 

nected in block djagonal matrix whose logical 
operators are {X^,- ■ ■ ,Xl^,Zl,- ■ ■ ,Zl^}, 
{Jq'^••• The blockcode is a 

n — ki code whose weight of elements should be 

redefined: as long as there has at least one ax or or ay 
in the first ki qubits we count " 1" , and the same to next 
k2, Kz, ■ ■ ■ , km qubits. We nest blockcode and subcode 
as the following: any physical qubit of blockcode is re- 
placed by the corresponding logical operator of subcode, 
then if we could find a redefined ^[Z^^iki,k' ,d!] code 
and if d is defined as the minimum sum of d' numbers of 
{di, ■ • • , dm}, a new [X^IILi k' , > d] degenerate QECC 
is constructed. 

The details will be presented in appendix. Our theory 
are general and for some special cases, it reduces to some 
known methods of constructing codes: (1) if all subcodes 
are the same codes and fc = 1, it will turn to the theory 
of concatenated coding (see IB2I) : (2) if the number of 
subcode is 2, it will turn to the theory of general stabilizer 
pasting for all distance . 
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APPENDIX A: CODES CONSTRUCTION FOR DISTANCE 3 



An [n, n — k,3] stabilizer code has k generators which build a (fc x 2n) generator matrix. There are 3n one-qubit 
error syndromes in which any bit flip error syndrome /(o",j) is the ith column of the right-half matrix, and any sigh 
flip error syndrome /(a* ) is the ith column of the left-half matrix and /(cr^) is the sum of these two columns. If the 
code is nondegenerate, then all the 3n syndromes are different. 

For a [n,n — k,d\ code to correct t errors, we define the using rate of syndromes: 

g{n,k,t)= ^^='J^'\ (Al) 

which is the division of the syndromes that all errors used and the syndromes that the code can provide. Using rate 
of syndromes is a very important parameter in codes construction, and obviously when g{n, k, t) < 1, the larger g is, 
the more optimal the code is. For a distance 3 nondegenerate code if g = 1 — we call it perfect code. 

There are four steps to construct a [n, n — k, 3] nondegenerate code: (i) choose 2n different syndromes for all Ux 
and az errors; (ii) calculate ay syndromes and make sure that all 3n syndromes are different; (iii) build the generator 
matrix with ax and a^ error syndromes and make sure that all k generators are commute to each other (we all it 
commuting condition): (iv) make sure that all generators are real generators which means no one is product of other 
generators (we call it gciuerator condition). Here we do not only restrict our focus on real stabilizer code but also 
discuss the complex stabilizer code, which means the number of ay in each generator need not to be even only. 

It is not easy to construct a code since step (iii) is a strong constraint condition and not easy to achieve. So 
constructing an optimal or near optimal code whose using rate of syndromes g{n, k, t) is equal to or near one is much 
harder. 



1. Nesting two [5,1,3] codes together 



[5, 1, 3] is the first prefect distance 3 code and has many good features that we need: (i) nondegenerate, which means 
all error syndromes are different; (ii) perfect, 5(5,4, 1) = 15/16, which means it used all the syndromes that the code 
can provide except for 0000; (iii) commuting, which means when we use [5, 1,3] codes to construct new codes, some 
parts of the commuting condition of new codes are satisfied automatically. 

A [25, 17, 3] code has 8 generators, which means the syndrome is a 8-bit number whose first 4-bit could define 
2^ = 16 syndromes. It is quite obvious that we can choose a [5, 1,3] code to take the first 4-bit syndromes, and the 
using rate will be 15/16. Here we call this [5, 1, 3] code blockcode. So the blockcode which takes the first 4-bit of the 
8-bit syndromes defines 15 blocks, i.e. {B^., B^, B^, B^, B^, B^, B^, B^, Bj, Bl,Bl, Bl,B'^, B^, B"^}, which are called 
blockcode syndromes and satisfy: 



Bl + Bl 



Bl 



(A2) 



and make all blocks different from each other. 

Then what we need to do is only to make the last 4-bit syndromes in any block different and make sure that any 

ay syndrome is the sum of corresponding ax and cr^ syndromes, then step (i) and (ii) of code construction will be 
finished. We can see that also a [5, 1, 3] code could satisfy this condition. We call this code subcode and the last 4-bit 



syndromes subcode syndromes, which is {S^, S^, S^, S^, S^, S^, S^, S^, S^, S^, Sl,Sy, Sy, Sy, Sy} and 

Q'i I Q'i Q'i 

So blockcode and subcode are nested in this way: 

{{Bl, Bl, Bl, Bl, Bl) ® {Sl, Sl, Sl, Sl, Sl), 

{bI,bI,bI, b1, bI) ® {Sl, sl, sl, sl, sl), 

{^l^Bl^By^^t'^y^ ® {'S'y,'S'y,5'^,S'^,5'^}}, 
which means if the [5, 1,3] blockcode and subcode are 



(A3) 



/ ail ai2 ai3 (ii4 ai5 

021 <J22 (l23 (l24 (l25 

0.31 032 033 034 035 

041 042 043 044 045 



^'11 ^12 ^'13 bi4 bi5 \ 

621 622 &23 ^24 b25 

&31 ^32 ^33 ^34 &35 

641 642 643 644 &45 / 



(A4) 



(A5) 



5 



the new nested code will be 



(Cll Ci2 Ci3 Ci4 Ci5 

C21 C22 C23 C24 C25 

C3I C32 C33 C34 C35 

C41 C42 C43 C44 C45 





di2 


du 


di4 


di5 


•^21 




d23 


d24 


d25 


C^31 




d33 


d34 


d35 


d4i 


di2 


d43 


da 


d45 



(A6) 



/ ail ail ail an an ai2 ai2 ai2 ai2 ai2 ai3 ais ais ais ais ai4 ai4 ai4 ai4 ai4 ai5 ai5 ai5 ais ais 

a2i a2i a2i a2i a2i a22 a22 a22 022 022 023 023 a23 a23 a23 a24 a24 a24 a24 a24 025 025 025 a26 a26 

a3i a3i a3i a3i a3i a32 a32 a32 a32 a32 033 ass ass ass 033 034 034 as4 as4 as4 035 035 035 ass ass 

a4i a4i a4i 041 041 a42 a42 a42 042 042 043 a4s a4s a4s 043 044 044 a44 a44 a44 045 045 045 a4s a4s 

Cll C12 CIS Cl4 Cl5 Cll C12 CIS Cl4 ClB Cn C12 CIS Cl4 Cls Cn C12 Cl3 Cl4 CIS Cll C12 Cl3 Cl4 Cls 

C2I C22 C2S C24 C25 C21 C22 C2S C24 C2B C21 C22 C2S C24 C2S C21 C22 C23 C24 C2S C21 C22 C23 C24 C2S 

CSl CS2 CSS C34 C35 C31 CS2 CSS CS4 C36 C31 C32 CSS CS4 CSS C31 C32 C33 CS4 CSS CSI C32 C33 C34 CSS 

\ C41 C42 C4S C44 C4B C41 C42 C4S C44 C4B C41 C42 C4S C44 C4S C41 C42 C43 C44 C4S C41 C42 C43 C44 C4S , « „s 

&11 611 bii 611 611 bi2 612 bi2 612 bi2 &1S 613 &1S 613 &1S 614 bi4 bi4 614 bi4 615 &is 615 615 &is \ 

&21 621 ^21 621 ^21 &22 622 &22 ^22 &22 f'23 f>23 &23 f>23 &23 f>24 &24 &24 ^24 ^24 ^25 ^25 f>25 625 &2S 

631 ^31 631 631 631 632 632 &32 ^32 ^32 ^33 hs hs hs ku ku ku ku ku kir, kib kir, &35 ^35 

641 641 641 641 641 642 642 642 642 642 643 &43 biS 643 643 644 644 644 644 644 645 645 645 &45 645 

(111 (ii2 dm di4 (ii5 dii di2 (ii3 ^14 dir, du drz din di4 dir, du dyz dm di4 dir, du dyz dyj, dn dm 

d21 d22 d23 d24 ^25 £^21 c!22 £^23 '^24 £^25 d21 ^22 £^23 £^24 £^25 £^21 £^22 £^23 £^24 £^25 c!21 £^22 £^23 £^24 £^25 
£^31 £^32 £^33 £^34 £^35 £^31 £^32 £^33 £^34 £^35 £^31 £^32 £^33 £^34 £^35 £^31 £^32 £^33 £^34 £^35 d3i d:i2 £^33 £^34 £^35 
£^41 £^42 £^43 £^44 £^45 £^41 £^42 £^43 £^44 £^45 £^41 £^42 £^43 £^44 £^45 £^41 £^42 £^43 C?44 £^45 £^41 £^42 £^43 £^44 £^45 / 

Step (iii) and (iv) are usually the hardest part for code construction, but in our method it is much easier to achieve. 
No matter how many times blockcode and subcode repeat, commuting conditions of themselves are also satisfied, so 
what we need to do is to satisfy commuting condition between blockcode and subcode which is 

{an + ai2 + ai3 + an + ai5){dji + dj2 + dj3 + dji + dj^) 
+{bii + bi2 + bis + bi4 + bi5){cji + Cj2 + cjs + Cj4 + Cjs) = (A8) 

then we can consider a more stronger constraint condition which is (a,i + a,2 + Oi3 + ai4 + Ois) =0 and (cji + Cj2 + 
•^js + ^j4: + ^35) = that means the number of "l"s in any row of left halves of generator matrices of both blockcode 

and subcode is even. 

For step (iv), obviou.sly blockcode and subcode satisfy generator condition themselves. It is easy to prove that 
between blockcode and subcode the only opportunity that one generator is product of others is (00 • • • 0|11 • • • 1) or 
(11 • • • 1|00 • • • 0) or (11 • • • 1|11 • ■ • 1) belongs to stabilizer of both codes at the same time, and the probability is quite 
small. For computer search Then we tiun to find codes that all (00 • • • 0|11 • • ■ 1), (11 • • • 1|00 • ■ • 0) and (11 • ■ • 1|11 ■ • ■ 1) 
do not belong to S which is a more stronger constraint condition for computer search. 

Depending on the computer search, we find many accordant [5,1,3] code, one of which is presented in the article, 
then the [25, 17, 3] new code by nesting two [5, 1, 3] codes is 

/XXXXXI I I I IZZZZZZZZZZYYYYYX 
ZZZZZXXXXXI I I I IZZZZZXXXXX 
ZZZZZZZZZZYYYYYI I I I lYYYYY 
I I I I IZZZZZZZZZZYYYYYXXXXX 
XIZZYXIZZYXIZZYXIZZYXIZZY ' 
ZXI ZXZXI ZXZXI ZXZXI ZXZXI ZX 
ZZYIYZZYIYZZYIYZZYIYZZYIY 

VlZIYXIZIYXIZIYXIZIYXIZIYX/ 

Here we discuss the using rate of syndromes in aspect of code nesting, g of both blockcode and subcode are 15/16, 
but each blockcode syndrome only meets 1/3 of all subcode syndromes. So the using rate is: 

1 75 1 

5(25, 8, 1) = -guocki^, 4, l)gsub{5, 4, 1) = ^ < 2 (^10) 

which is the same as the result that calculated by equation (1). As we know that the optimal 25-qubit code is [25, 18, 3] 
whose using rate is 75/128, this code is near optimal. So g is an important parameter to judge whether a code is 
optimal or not: generally if | < 3 < 1, the code must be optimal. 

In the same way, we nest this [5, 1,3] code and this [25, 17,3] code to construct a [125, 113,3] code, and it is easy 
to prove that the code satisfies all constructing conditions automatically. So all the [5", 5" — 4n, 3] codes could be 
constructed. 
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Then the [30, 22, 3] near optimal code takes the form: 

{{Bl, BI, BI, B^, 5^,0000} ® {sl, sl, s^, sl, sl), 
{Bl,Bl, Bl, Bf, Bl, 0000} ® {Sl,Sl, Sj, S'f }, 
{Bl,Bl, Bl, Bl, Bl, 0000} ® {S^ 5^, 5*, Sl}}, 

so the generator matrix is 

/XXXXXI I I I I ZZZZZZZZZZYYYYYI I I I I\ 

ZZZZZXXXXXI I I I IZZZZZXXXXXI I I I I 

ZZZZZZZZZZYYYYYI I I I lYYYYYI I I I I 

I I I I IZZZZZZZZZZYYYYYXXXXXI I I I I 

ZYXI ZZYXI ZZYXI ZZYXI ZZYXI Z 



X I Z 

ZXI ZXZXI ZXZXI ZXZXI ZXZXI ZXZXI 
ZZYIYZZYIYZZYIYZZYIYZZYIYZZY 
lYXIZIYXIZIYXIZIYXIZ 



Y 
X 
Y 



\ I Z I Y X I Z 



I Y X/ 



Similarly all [6 • 5", 6 • 5" - 4(n + 1), 3] code and [5 • 6", 5 • 6" — 4(n + 1), 3] codes could be constructed. 
Then the [35, 27, 3] code from the [30, 22, 3] code takes the form: 

{{^1, Bl, Bl, Bl, Bl, 0000} ® {Sl, sl, Sl, Sl, Sl}, 

{Bl,Bl,Bl,Bl,Bl}®mm, 

{Bl, Bl, Bl, Bl, Bl 0000} ^ {Sl,Sl, Sl, Sl, Sl}, 

{Bl,BlBl,Bl,Bl}®mm, 

{B]/, B^, Bl, B^, B^, 0000} (g) {Sl,Sl, Sl, S^, S"^}, 
{Bl,Bl,Bl,Bl,Bl}(i^ 0000} 

whose generator matrix is 



(All) 



(A12) 



(A13) 
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(A14) 



2. Constructing optimal codes 

In our method to construct a distance 3 code, since each blockcode syndrome only meets 1/3 of all subcode 

syndromes, the using rate 

9 = ^guockgsub (A15) 

If a code is nested by two optimal codes, g will always be less than 1/3, which means the code will never be optimal. 
So if we want to construct an optimal code, we need to find codes whose g is more than one to be blockcode or 
subcode. We call these codes whose g are more than one over — optimal codes. 
Suppose block and subcode are nested in this form: 

{{Bl-- - ,B:i}®{si,--- ,s:'}, 
{Bl,.-. ,b:}®{si,--- ,s:'}, 

{Bl,... ,B^}®{Sl,... ,S;;'}}, (A16) 

and suppose blockcode is an optimal distance 3 code, since each B^ syndrome only meets all S^ syndromes and each 
Bz syndrome only meets all Sz syndromes and each By syndrome only meets all Sy syndromes, there is no need to 
make any syndrome of subcode be different, which means the distance od subcode need not to be 3. We just need to 
find out a code whose Sx, Sz and Sy syndromes are different themselves. 
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a. the optimal [10, 4, 3] QECC 



As we know, the optimal 10 qubits distance 3 code is [10,4,3], which means if we want to construct this code, the 
subcode must be a [2,0, < 3] over-optimal code which we write as {2,2}. Suppose blockcode and subcode are 



/ Cn Ci2 Ci3 Ci4 Ci5 

C21 C22 C23 C24 C25 

C31 C32 C33 C34 C35 

\ C41 C42 C43 C44 C45 



dii di2 di3 di4 dis 

(^21 C^22 ^23 C^24 C?25 

(^31 <i32 C^33 C^34 ^35 

dil d42 d43 d44 d45 



(A17) 
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then the new nested code will be 
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Easily we can see that the commuting condition is 

(cjihix + c-j\bi2 + Cj2bt\ + Cj2bi2 + Cj36ji + Cj^ha + Cj^fe^ + Cj4&i2 + c^s^ii + 0^-5612) 
+((ijiaii + djiai2 + dj2aii + dj2ai2 + djsUn + dj3ai2 + dj4aii + dj4ai2 + dj^an + dj5ai2) 
= {cji + Cj2 + Cj3 + Cj4 + Cjr,){ba + &i2) + (dji + dj2 + djs + dj4 + dj5){aii + 0^2) = 0. 



(A18) 



(A19) 



(A20) 



We also consider more stronger constraint conditions instead, and different constraint conditions will give different 
constructions. Actually we should not give over-optimal code too many constraint conditions since the code may be 
quite difficult to be found or even do not exist. So we consider the conditions as (ci + C2 + C3 + C4 + C5) = and 
{di +d2 + d3 + d4 + ds) = which are both for the [5, 1, 3] blockcode. Then by computer search we have the blockcode: 



(A21) 
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and the subcode: 
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(A22) 



So the generator matrix of this [10, 4, 3] code is 
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(A23) 



6. the optimal [15, 9, 3], [20, 13, 3] and [25, 18, 3] QECC 



Similarly we can find {3, 2}, {4, 3} and {5, 3} over-optimal codes which are 
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(A24) 
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so the new [15, 9, 3], [20, 13, 3] and [25, 18, 3] codes are: 
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(A27) 





X 


X 


X 


I 


I 


I 


I 


z 


z 


z 


z 


Z 


z 


z 


z 


Y 


Y 


Y 




z 


z 


z 


z 


X 


X 


X 


X 


I 


I 


I 


I 


z 


z 


z 


z 


X 


X 


X 


X 


z 


z 


z 


z 


z 


z 


z 


z 


Y 


Y 


Y 


Y 


I 


I 


I 


I 


Y 


Y 


Y 


Y 


I 


I 


I 


I 


z 


z 


z 


z 


z 


z 


Z 


z 


Y 


Y 


Y 


Y 


X 


X 


X 


X 


X 


I 


z 


X 


X 


I 


z 


X 


X 


I 


Z 


X 


X 


I 


Z 


X 


X 


I 


z 


X 


z 


X 


z 


Y 


z 


X 


z 


Y 


z 


X 


z 


Y 


z 


X 


z 


Y 


z 


X 


z 


Y 


\I 


z 


X 


Y 


I 


z 


X 


Y 


I 


z 


X 


Y 


I 


z 


X 


Y 


I 


z 


X 


yJ 



(A28) 
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c. Another construction of [10, 4, 3] code 

Here we give another way of constructing [10,4,3] code. We may notice that two copies of a [5,1,3] code whose 
syndromes are {{/S^}, {S'^), {S'y}} could build a {10, 4} over-optimal code in this way: 

{{S.} = [S',} u [S',}, 
{S.} = {S'^}^{S'y}, 

{Sy} = {S'y} U {5;}}, (A30) 



then wo choose 



as blockcode, so commuting conditions of blockcode and subcode are satisfied automatically. 



Commuting conditions between blockcode and subcode is ai(di + • • • + dio) + &i(ci + • • • + cio) = 0, which means 
the number of "l"s in any row of both halves of generator matrices of {10,4} code is even. Instead we turn to 
search a [5, 1, 3] code which satisfies the same condition, then it is easy to prove that the constructed {10,4} code is 
appropriate. It is obvious that the [5, 1,3] code above 



/ 1 1 

10 1 

10 10 

V 1 1 



1 1 0\ 

110 

1 1 

1 l/ 



(A31) 
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is feat. Then the [10, 4, 3] code is 
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(A32) 



3. Gottesman's [2*, 2*^ — fe — 2, 3] codes and perfect codes 

The over-optimal codes whose Sx, Sz and Sy syndromes are different themselves are very important codes in our 
method of constructing optimal codes, so we call them opfAmal — constructing codes. More important thing is to 
know g's upper and lower bounds of optimal-constructing codes and when the codes could reach the upper bound. 
Here we have the lemma below: 

Lemma: For a [n,n — k, d\ optimal-constructing code, I < g < — - < 3. When n = 2*^ — 1, the code could reach 
the upper bound. 

Proof: A [n.n — k,d] code could provide 2^^ syndromes. If each of Sx, Sz and Sy takes all the syndromes except 

for 00 •••0 and S^ + S^ = S^, this code must reach the upper bound. So g's upper bound is — '-. But not 
optimal-constructing codes of all length could reach upper bound, since each of Sx, Sz and Sy takes all the syndromes 
except for 00 • • • 0, n must equal to 2*^ — 1. 

So [2'' — 1,2^^ — 1 — fc,< 3] codes are optimal-constructing codes that reach the upper bound, and we call them 

raw — perfect — constructing codes. Since limfc^oo - 2^^ — ^ ^' these code are quite useful to constructing optimal 
codes and perfect codes. And for the raw perfect-constructing codes wc have the lemma below: 

Lemma: Suppose any generator of a [2*^ — 1, 2*^ — 1 — /c, < 3] code is (ci, • • • , C2k_i\di, • • • , rf2'=-i), then (ci H h 

C2fc_i) = and (di -|- • • • -|- d2k_i) = 

Proof: Since this code used all syndromes that the code could provide, any row of the left half of generator matrix 

will have 2'°"^ "l"s which means the number of "l"s is always even, so (ci H h C2k_i) = 0. The same to the right 

half which means {di + h ^2*=-!) = 0. 

1 . Gottesman's [2'^,2'^ — fc — 2,3] codes are one of this kind. A [2^ , 2^^ — fc — 2, 3] code of this class can be constructed by 

nesting a "code" which is ^ q ^ and a [2^^ — 1, 2*^ — 1 — fc, < 3] raw perfect-constructing code together. ^ ^ ^ 

is blockcode whose syndromes are {10,01,11} and the raw perfect-constructing code is subcode whose syndromes 
a.re:{{Sx}, {Sz}, {Sy}}- The code is constructed this form: 

{{10}®{Sx,00---0}, 
{01}®{5„00---0}, 

{11} (g) {5^,00 •••0}} (A33) 

So g of this class is always 5 ■ | • 3 = | . 

Here we show how to construct the second code of this class- [16, 10,3]. Firstly we need to construct {15,4} raw 
perfect-constructing code. We could choose any [5, 1, 3] code whose syndromes are {{-S*^}, {S'z}, {Sy}}, and the {15, 4} 
code takes this form: 

{{Sx} = {S'x} U {S'z} U {S'y}, 
{Sz} = {S'z} U {S'y} U {S'x}, 

{Sy} = {S'^}U{S'JU{S'J} (A34) 

It is easy to prove that SI.+SI = Sy and the commute conditions of subcode and blockcode are satisfied automatically. 

Commute condition between subcode and blockcode is: ai(di -|- h die) + &i(ci H h cie) = 0, and according to 

the lemma above this condition is also satisfied automatically. So if we choose a [5, 1, 3] code as 



/X I Z Z Y\ 

Z X I Z X 

Z Z Y I Y 

\I Z Z Y xj 



(A35) 
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then the [16, 10, 3] code is 
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(A36) 



,3,3] code have been constructed which is 



1 1 1 1 1 1 1 1 
00000000 
1 1 1 1 
1 1 1 1 
VO 1 1 1 1 



0000000 0\ 
11111111* 
10 110 10 
10 10 10 1 
110 10 1/ 



(A37) 



But here we give another construction of [8, 3, 3] that out of the method of building [2^^, 2*^ — fc — 2, 3] codes, which is: 



/ 1 1 1 1 

1 1 1 1 

1 1 1 1 

1 1 1 1 

\0 1 1 1 1 



1 1 1 1 0\ 
lOlOOlOl' 
110 10 
10 10 1 
10 1/ 






1 1 



1 
1 



(A38) 



2. The perfect one-error-correcting codes are also one of this kind. This class is [(2*^ — l)/3, (2'^' — l)/3 — fc, 3] codes 



whose k must be even. Though the first one of this class which is 



1 
1 



is not a real code, we still call it " perfect 



code" since it is quite important in code construction. Comparing with other methods of constructing perfect codes, 
our method is more general. We have the theory below: 

Theory: A [(2'=+'=' - l)/3, (2'=+'=' - l)/3 - (fc + fc'), 3] perfect code can be constructed by nesting a [(2'= - l)/3, (2'= - 
l)/3-/c,3] perfect code, a [2^=' - 1, 2'^' - 1 - fc', < 3] raw perfect-constructing code and a [(2*^' -l)/3, (2*^' - l)/3 - fc', 3] 
perfect code together. 

Proof: Firstly we nesting a [(2*^ - l)/3, (2'= - l)/3 - fc,3] perfect code and a [2'^' - 1,2'^' - 1 - fc', < 3] raw 
perfect-constructing code together in this way: 



{{B,}(S>{S^,00---0}, 
{SJ(»{5„00---0}, 
{By} ^{Sy, 00 ■■■0}} 



(A39) 



Since all commute conditions are satisfied automatically, this construct a [(2*^+*'' — 2''" )/3, (2'^'+'^ —2'^ )/3 — (fc -I- fc'), 3] 
optimal code. Notice that we never used syndrome 00 • • • of blockcode, and the raw perfect-constructing code 
has fc' generators, so we could use syndrome 00 • • • to nest another perfect code which also has fc' generators. 
This [(2'=' - l)/3,(2'=' - l)/3 - fc',3] perfect code is another subcode with syndromes {S'^, S'^, S'y}. Finally this 
pfc+fc' „ ly^^ (^2''+''' - l)/3 - (fc + fc'), 3] code is constructed: 



{{B,} {S^, 00 • • • 0}, {B,} {S,, 00 • • • 0}, 
{By} 8) 00 • • • 0}, 00 • • • ® {S'J, 
00---0(g>{S',},00---0(g>{S'}} 



(A40) 



It is easy too prove that all commute conditions are satisfied. 

In process of code construction, g of the [(2*^+'^ - 2'^')/3, (2*^+*^ 
i • (1 - ^) - 3 = 1 - ^ and5 of 00 ■•• pasted with [(2*=' - l)/3, (2'--' - 
the [(2'^+'=^' - l)/3, (2'=+'=' - l)/3 - (fc + fc'), 3] code g = gi + g2 = 1 - 

Actually our method is not constrained by whether fc and fc' are even or odd, which is mentioned in Ref.fo']. For 
example, [40,33,3] code can also be constructed using our method. 



2*^ )/3 — (fc -I- fc'), 3] optimal code is gi = 
l)/3 - fc', 3] code is 32 = ^ ■ (1 - ^), so for 
^ Obviously it is a perfect code. 
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4. Gottesman's stabilizer pasting 



Stabilizierpasting: Given two nondegenerate stabilizer codes [112, S2] — {Si,-- - jS'^j) and [711, si] = (Ti,--- ,Ts-^) 
of distance 3, if two observables X(n2) and Z{n2) belong to [rt2,S2], say, 5*1 ~ X{n2) and S2 = ^("■2), then the 
stabilizer constructed in the way below defines a nondegenerate stabilizer code [712 + ni, s] with s — max{s2, si + 2}, 
denoted as [712,52] > [7^l,Sl]. 



Consider his theory with our method, the blockcode of [712,32] is always ^ j wc mentioned above, and the 

subcode may not be a real code since commute conditions may not be satisfied. Notice that syndrome 00 of blockcode 
has never been used, so we put another code [711, si] under 00 to finish the construction. Code [771, si] may also not 
be a real code as long as the final constructed code satisfies all commute conditions. Many optimal codes have been 
constructed in Ref. pl| . 

To construct the [37, 30, 3] code, firstly we choose the nontrivial [8, 3, 3] code found above 
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(A41) 



as blockcode and a {4, 2} code found by computer search 
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(A42) 



to be subcode. There have 8 syndromes that the [8, 3, 3] code did not used which are: 



/O 1 1 1 1\ 
01010101' 
10 10 10 1 
110 110 

Vo 1 1 1 1 o/ 



(A43) 



Since syndromes of [37, 30, 3] code are 7-bit numbers, with another two bits we can find a {5, 7} code below: 








10 10 
10 10 
10 10 



10 10 
10 1 
10 10 

Vo 1 1 



1 1 \ 

1 1 

1 1 

110 

110 

1 1 

10001/ 



(A44) 



The last four generators build a [5, 1, 3] code and the syndromes of first five generators all belong to those 8 syndromes 
mentioned above. It is easy to check that all commute conditions are satisfied, then the [37, 30, 3] code is constructed 
as 



/X XXXYYYYXXXX 
ZZZZIIIIZZZZ 
IIIIZZZZIIII 
I I I IZZZZXXXX 
I I I lYYYYZZZZ 
YZIZYZIZYZIZ 

\Y X Z I Y X Z I Y X Z I 



YYYYZZZZI I I I 
I I I IXXXXYYYY 
ZZZZYYYYXXXX 
Y Y Y Y I I I I Z Z Z Z 
XXXXZZZZXXXX 
YZIZYZIZYZIZ 
YXZIYXZIYXZI 



Z Z Z Z I I I I Z X I X Z \ 
XXXXYYYYZXIXZ 
YYYYXXXXZXIXZ 
XXXXYYYYXZZX I 
I I I lYYYYIXZZX 
YZIZYZIZXIXZZ 
YXZIYXZIZXIXzJ 

(A45) 
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APPENDIX B: CODES CONSTRUCTION FOR ALL DISTANCE 



Suppose there are two copies of one [5,1,3] code whose generators and logical operators are Mi, M2, M3, M4 and X, Z 
pasted together. If the construction is ^ ^ , then the logical operators of this code are {X^I, I(^X, Z®I, I®Z}; 

if the construction is [A A ) , then the logical operators of this code are {X (^I,I^X,Z(S)I,I(S)Z, Mi (g) /, M2 ^ 
I, M3 (g) J, M4 (g) /, Ml (g) Ml, M2 (X" M2, (g) M3, M4 (g) M4}. Obviously the first one is the easiest way to construct 
good blockcode and the second one is the hardest. 



1. Concatenated coding 



Suppose 5 copies of a [5, 1, 3] code are pasted in the easiest way, then the constructed [25, 20] code has 20 generators 
and 10 logical operators. So if we want to build code of larger distance, we must pick some logical operators and put 
them into the stabilizer, and these logical operators build the blockcode. Suppose the logical X and Z operators of the 
five copies are Xi , X2 , X^ , X4 , X^ and Zi, Z2, Z^, Z^, Zq, it is easy to think out that a 5-qubit code whose physical 
qubits are replace by these ten logical operators can be constructed. We choose this 5-qubit code to be blockcode, 
and if distance of blockcode is d then distance of this constructed code will be 3d. As we know the largest distance of 
5-qubit codes is 3, so choose any [5, 1, 3] code to be blockcode we will get a [25, 24, 9] code. For example, if we choose 

/lOOlO OllOOX 

1 1 1 1 I , . 

10100 00011' 

VoiOlO 10001/ 

then the generator matrix is 
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(B2) 



This is so-called concatenated coding which has already been pointed by Gottesman in Ref. [j]. The codes con- 
structed in this way must be degenerate code since distance only increases in N{S) — S but distance of S does not 
change. Here we see that the concatenated coding construction can be understood in the framework of the nested 
coding construction. 
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2. Stabilizer pasting for all distance 

Theory of general stabilizer pasting for any distance is: 

Theory: Suppose there have four stabilizers Ri, i?2, 5*1, and ^2 with Ri C 5*1 and i?2 C 5*2 • Let Ri define a 
[ni,^i,Ci] code, i?i be a [n2,^2)C2] code, 5*1 be a [ni,ki,di] code and ^2 be a [n2,k2,d2] code with ki < li, Ci < di, 
h — ki = I2 — k2 and Si, S2 to be nondegenerate. Let generators of -Ri be {Mi, • • • ,M„^_;j}, generators of Si be 
{Ml, • • • , Mni-kA, generators of R2 be {A^i, • • • , iV„2_/2}, and generators of ^2 be {Ni, ■ ■ ■ , N^^-kA- Then we can 
build a new stabilizer S' on ni + n2 qubits generated by 

{Ml®/,--- ,M^,_i,®I,I®Ni,--- ,I^N„^_i,, 

M„,_,,+i ®7V„,_,2+i>-- - ,M„,_fc, ® iV„2-feJ- (B3) 

The code has (ni — h) + {n2 — I2) + {k — ki) generators which means it encodes ^1 + ^2 = ^2 + ki qubits, and the 
distance of the new code is min{di,d2, ci +02}- 

Since Ri and R2 are connected in block diagonal matrix, obviously they are subcodes and the blockcode is built 
by all logical operators ofthem. Suppose logical operators of Ri are {Xi,--- ,Xi^,Zi, - ■ ■ ,ZiA, and logical opera- 
tors of R2 are {X[,--- ,X[ ,Z[,--- ,Z[ }. Then we can rewrite Si as {Mi,--- ,M„j_;j,Xi, ■ - - and ^2 as 

{TVi,--- ,7V„2_;2,X(, - - - so logical operators of Si and S2 are {Xi^^kt+i, - ■ ■ , Xi^, Zi^^ki+i, ■ ■ ■ , Zi^ and 

{X'i^_^^^^,- ■ ■ ,X[^,Zl^_^^^^, - ■ ■ ,ZIJ. Obviously the blockcode is {XiiE)X[,--- ,Xi^^ki ® X[^_^A and the logical 
operators of the constructed code are {Xi^^k^+i (gi /, - - - , Xi^ ® I, Zi-^^ki+i <8> -f, • ' • , Zi^ 1^ 1} Li {I ® -'^;'2-fe2+i' " ' ' ® 

,i^z\}u{Xi(3i,--- ,Xl;_k^^I}^{Zl®z{,■■■ ,zi;_k,®z'~_^A. 

We can see that there are three kinds of logical operators. The first one are logical operators of ^i and S2, which 
act on stabilizer and make the minimal weight be min{di, ^2}- The second one are {Xi ® /, - - - , Xi-^^-ki ® and 
since Si and 6*2 are nondegenerate, they act on stabilizer and make the minimal weight be min{di, ^2}- The last one 
are {Zi ® Z[, ■ ■ ■ , Zi^^ki ® Z[ }, which act on stabilizer and make the minimal weight be Ci + C2. So the distance 
of the new constructed code is min{di, ^2, ci + C2}. 



